Non-Linear Effects in Resonant Tunneling; Bistabilities and Self-Sustained Oscillating 

Currents 
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We study non-linear phenomena in double barrier het- 
erostructures. Systems in 3D under the effect of an external 
magnetic field along the current and ID systems are analyzed. 
Non-linearities are reflected in the I-V characteristic curve as 
bistabilities, instabilities and time dependent oscillations of 
the currents. The nature of the non-linear behavior depends 
upon the parameters that define the system. 



I. INTRODUCTION 

Resonant tunneling in double barrier heterostructures 
(DBH) have been found to exhibit a variety of interesting 
physical propertiesH. Transport in a 3D DBH presents 
non-linear phenomena, which are reflected in the obser- 
vation of multistabilitics in the I-V characteristic curvea. 
This property is produced by the rapid leakage of the 
electronic charge accumulated at the well between the 
barriers when the applied potential is just taking the de- 
vice out of resonance. It is a result of the Coulomb charge 
interaction at the well and has been theoretically stud- 
ied, in a 3D system, using a mean field approximationp. 
In this work we have found a similar behavior in a ID 
system constituted by a quantum dot connected to two 
leads. Although, the origin of the bistability is the non- 
linearity produced by the Coulomb interaction, as in the 
3D system, what is dominant in a quantum dot is the 
spatial e-e correlation. It produces simultaneously the 
Coulomb blockade of the charges, that appears as peri- 
odic oscillations of the conductance as a function of the 
gate potential and the bistability of the current. For a 3D 
system instead, the Coulomb blockade is absent and the 
electronic interaction can be incorporated as a renormal- 
ization of the quasi-particle energies using a mean field 
approximation which gives rise to the non-linear effects. 

Non-linear time dependent transport in double-barrier 
structures is another very interesting phenomena. In this 
work we report the presence of self-sustained current os- 



cillations when a large magnetic field is applied along the 
direction of the current in a 3D DBH. The magnetic field 
produces a drastic change in the number of resonant elec- 
trons going through the well inducing a great enhance- 
ment of the non linearities of the system. The state of 
charge of the well controls the oscillatory behavior of the 
current. Our results show that the system bifurcates as 
the magnetic field is increased and may transit to chaos 
at large enough field. The oscillations are exponentially 
sensible to initial conditions or to infinitesimal changes 
in the parameter of the Hamiltonian, characterizing a 
chaotic behavior which has been studied in detail in the 
parameter space. 

We obtain a similar result for the current going along 
the quantum dot connected to two leads. The produced 
self-sustained time dependent current oscillations is con- 
trolled in this case by the voltage applied to the metallic 
gate. In both cases the period of the oscillations depends 
upon the transit time of the electron going along the well, 
which is determined by the size of the barriers. 



II. THE MODEL 



Consider a ID or a 3D double barrier structure (DBS) 
device. In order to study its time evolution under a bias 
we adopt a first-neighbors tight-binding model for the 
Hamiltonian. For the 3D situation we include an exter- 
nal magnetic field B in the growth direction. We assume 
a parabolic energy dispersion parallel to the interfaces 
while the field quantizes the motion of the electrons in 
the perpendicular plane giving rise to Landau levels with 
energies e„ = (n + h)hcj c , where n = 0, 1, 2, ... is the Lan- 
dau index and u c = eB/m*c is the cyclotron frequency. 
The probability amplitude bj for an electron in a time 
dependent state \kn > at plane j along z in the Landau 
level n is determined by the equation of motion 
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where v is the hopping matrix elements between electrons 
in nearest-neighbor planes, and ej represents the band 
contour and external bias. Here the sum over (m, k) cov- 
ers all occupied electron states of the system, within the 
various Landau levels m with energies below the Fermi 
level, incident on the DBS from the emitter side. In 
writing equation [??] we have adopted a Hartree model 
for the electron-electron interaction. We have neglected 
inter-Landau level transitions since they_are of little im- 
portance at not too low magnetic fieldstj, and averaged 
over allowed transitions between degenerate states within 
a Landau level, taking advantage of the localization in- 
duced by the Gaussian factor in the Landau basis. 

For the ID system we can write a similar equation 
where the subindex n corresponding to the Landau level 
is eliminated. The electron-electron interaction restricted 
to the dot has to be treated in this case using a strong in- 
teraction limit approximation capable of including simul- 
taneously the Coulomb blockade effect and the non-linear 
behavior. The electron-electron interaction is analyzed 
within the context of the Hubbard I approximation! In 
this case the equation of motion is written as, 

db k 

fo-jr = + + (» ¥= 0); 



in aT = y° + Vp ) b i a + nav ° y- 1 + 5 v (2) 

where a = ±; n + = n; n~ = 1 — n; e + = U; e~ = 
0; v 3 =v Vj £ 0,1; Vj = v = 0, 1; n = £ fe , \b k '\ 2 ; 

b k = J2a=i bo a 'i Vp m tne gate potential applied at the 
well and the sum over k' covers as before all occupied 
electron states of the system. 

The time dependent equations are solved using a half- 
implicit numerical method which is second-order accu- 
rate and unitaryu. 

The approach taken here assumes that in the presence 
of a bias, the wave function at time t outside the structure 
is given byEI 

b nk = Ie} kr 3 + R . ne -ikr } f < _ £ 

bf = T jn e lk ' r > r > -L (3) 

Here k and k' = y / '2rn*[e k + e£\/Ti are the wave numbers 
of the incoming and outgoing states, respectively, with 
e k = e n — 4vsin 2 (^) the energy of the incoming parti- 
cle (for the ID case e„ = 0) and cl is the external bias. 
To model the interaction with the particle reservoir out- 
side the structure the incident amplitude I is assumed to 
be a constant independent of the coordinates. The en- 
velope function of the reflected and transmitted waves, 



Rj n and Tj n , are allowed to vary with j, however. Since 
far from the barriers these quantities are a weak func- 
tion of the coordinate zj we restrict ourselves to the lin- 
ear corrections only. This approximations is appropriate 
provided the time step St does not exceed a certain lim- 
iting value. For the results presented here, a maximum 
value of St = 3 x 10 _17 s was found sufficient to eliminate 
spurious reflections at the boundary while maintaining 
numerical stability up to 20 x 10~ 12 s. 

In our numerical procedure the coefficients obtained 
for one bias are used as the starting point for the next 
bias step. 

We study first the 3D DBH under the effect of an ex- 
ternal magnetic field. The model is applied to an asym- 
metric double barrier structure. The second barrier is 
made wider than the first in order to enhance the trap- 
ping charge in the well. The conduction band offset is 
set at 300meV. In equilibrium and at B = 0T, the Fermi 
level lies 19.2meV above the asymptotic conduction band 
edge, so that the zero bias resonance lies well above the 
Fermi sea. The parameter values in equation [??] are 
set at v = — 2.16cV and U = lOOmeV. The latter was 
chosen phenomenologically so as to fit the experimental 
I-V characteristic for apGaAs device in the absence of an 
external magnetic fielcfl. 

For small values of the magnetic field a stationary so- 
lution is reached after a short transient. At large enough 
bias two stationary solutions emerge, obtained by cither 
increasing or decreasing the applied voltage. This corre- 
sponds to the well known bistability of 3D DBH. 

A completely novel feature starts up as the field is in- 
creased. At small values of the external bias a stationary 
solution is rapidly reached. As the bias is increased, how- 
ever, a critical value arises beyond which no stationary 
solution exists, and the system begins to oscillate in a 
self-sustained way. After a range of voltages over which 
the oscillations persist a stationary solution is reached 
again. 

In fig.l we show the current at the center of the well 
for B = 13T and a fixed bias V = 0.27v. A transient 
period is followed by an oscillation that is never damped 
out. Although not perfectly periodic, the oscillation has 
two strong Fourier components at frequencies v ~ 0.3 
and 0.8THz. 

For a ID system we study a model which consists 
of leads connected to a quantum dot represented by 
v = O.lv, U = l.Ov and a Fermi level Ef = O.Olw. The 
normalization of the wave function is taken so that each 
site could have up to a maximum of two electrons. We 
solve first equations || supposing that they possess sta- 
tionary solutions. Since the electron density has rather 
long-range oscillations we made sure the sample was long 
enough to make finite size effects negligible. This is guar- 
anteed taken a sample of N > 100 sites. Fig 2a shows 
the current voltage (I-V) characteristics for several a gate 
potential applied at the quantum dot V p = O.lv. The 
system has a bistable behavior, which corresponds to the 
well known two solutions for the current; as the voltage is 
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increased when the dot is charged (solid line) , and when 
the voltage is decreased the solution when the dot is with- 
out charge (dashed line) . In cither case the self-consistent 
solution converges to a stable fixed point after some finite 
number of iterations. Similarly to 3D tunneling devices in 
this ID structure an increase of the Fermi level augments 
the charge content of the well and consequently the non- 
linearities and in particular the width of the hysteresis 
loop. When the gate potential is greater than a thresh- 
old value, which depends upon the parameters of the sys- 
tem, a completely novel feature starts to develop where 
two periods, two fixed points are encountered for certain 
voltages as shown in Fig 2b. This bubble-like solution 
is obtained as the voltage is increased outside the region 
of bistability showing that although both phenomena are 
derived from the non-linearities introduced by the local 
Coulomb repulsion, they corresponds to two different ef- 
fects. Note that these new solutions are not reached, as 
the voltage is decreased (dashed line) as it corresponds 
to a discharged dot where the non-linearities are negli- 
gible. As the potential gate is raised still further, the 
area of the bistable bubble increases and it enters into 
de bistable region. At larger gate potentials the solution 
undergoes further bifurcations and finally a bifurcation 
cascade leading to a chaotic region. These behavior is in 
fact indicating the absence of stationary solution of the 
system. 

Similarly to the 3D DBH under the effect of a mag- 
netic field, self-sustained oscillations for the current are 
found in the case of a dot connected to two leads in the 
parameter region where there is no stationary solutions. 

We interpret these results for both systems in the fol- 
lowing way. In general, current flows through the het- 
erostructure as long as a tunneling resonance lies within 
the emitter Fermi sea. We assume that at very low bias 
the resonance lies above the Fermi energy so that no 
current flows. As the bias increases it drops, reaching 
eventually the Fermi energy, thus opening a channel for 
electrons to tunnel through the double barrier. A current 
is thus established. The charge is trapped in the poten- 
tial well, raising its bottom and pushing the resonance 
towards higher energies for the 3D system. For the ID 
system the weight of the resonance level is proportional 
to 1 — n and as a consequence reduces when the charge 
increases. The current drops, some of the accumulated 
charge leaks out allowing the current to flow more easily 
once again, and a new cycle begins. 

According to the picture drawn above the charge in 
the well lags the current, as exhibited in the inset of 
fig.l. Here the charge at the center of the well (dashed 
line) is plotted together with the current at the same 
point (full line), the former displaced a time r ~ 1.4ps to 
the left with respect to the latter. The relaxation time r 
and the period of the oscillations are determined by the 
tunneling time for the electrons to leak out through the 
barriers and may be adjusted by modifying the barrier 
thickness. 



III. CONCLUSIONS 

There is a wide region of accessible values of these pa- 
rameters were bistability, instabilities and oscillations in 
the current are either present or absent. This opens up 
the interesting possibility of experimentally studying the 
transition from one behavior to the other as the param- 
eters of both systems studied are varied. From the view 
point of possible applications these systems could operate 
as electromagnetic generators in the THz regime. 
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FIG. 1. Current as a function of time at the center of the 
well for a 3D DBH under the effect of an applied magnetic 
field 13T. The inset shows the charge and the current shifted 
one relative to the other by 1.4ps. 
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FIG. 2. I-V 

characteristic of a quantum dot connected to ID leads for: 
a) V p = O.lv and b) V p = 0.2v 
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